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Abstract 

Higgs boson decays mediated by flavor changing neutral currents (FCNC) are 
very much suppressed in the Standard Model, at the level of 10 -15 for Higgs boson 
masses of a few hundred GeV. Therefore, any experimental vestige of them would 
immediately call for new physics. In this paper we consider the FCNC decays 
of Higgs bosons into a top quark in a general two-Higgs-doublet model (2HDM). 
The isolated top quark signature, unbalanced by any other heavy particle, should 
help to identify the potential FCNC events much more than any other final state. 
We compute the maximum branching ratios and the number of FCNC Higgs boson 
decay events at the LHC collider at CERN. The most favorable mode for production 
and subsequent FCNC decay is the lightest CP-even state in the Type II 2HDM, 
followed by the other CP-even state, if it is not very heavy, whereas the CP-odd 
mode can never be sufficiently enhanced. Our calculation shows that the branching 
ratios of the CP-even states may reach 10~ 5 , and that several hundred events could 
be collected in the highest luminosity runs of the LHC. We also point out some 
strategies to use these FCNC decays as a handle to discriminate between 2HDM 
and supersymmetric Higgs bosons. 
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At the tree-level there are no Flavor Changing Neutral Currents (FCNC) processes in the 
Standard Model (SM), and at one- loop they are induced by charged-current interactions, 
which are GIM-suppressed [I]. Letting aside the meson-meson oscillations, such as K° — 
K° and B° — B°, the decay processes mediated by FCNC are also of high interest and 
are strongly suppressed too. For instance, we have the radiative B-meson decays, with a 
typical branching ratio BR(b —>■ S7) ~ 10 -4 . But we also have the FCNC decays with the 
participation of the top quark as a physical field, which are by far the most suppressed 
decay modes j2MHl- Indeed, the top quark decays into gauge bosons (t ^ cV; V = 7, Z, g) 
are well known to be extremely rare events in the SM. The branching ratios are, according 
to Ref. [Sj: ~ 5 x 10~ 13 for the photon, slightly above 1 x 10~ 13 for the Z-boson, and 
~ 4 x 10~ n for the gluon channel, or even smaller according to other estimates j3]. 
Similarly, the top quark decay into the SM Higgs boson, H SM , is a very unusual decay, 
typically BR(t — > cH SM ) ~ 10~ 14 [5j. However, when considering physics beyond the 
SM, new horizons of possibilities open up which may radically change the pessimistic 
prospects for FCNC decays involving a Higgs boson and the top quark. For example, 
in Ref. jH] it was shown that the vector boson modes can be highly enhanced within 
the context of the Minimal Supersymmetric Standard Model (MSSM) [7j. This fact 
was also dealt with in great detail in Ref. [S] where in addition a dedicated study was 
presented of the FCNC top quark decays into the various Higgs bosons of the MSSM (see 
also 0), showing that these can be the most favored FCNC top quark decays - above the 
expectations on the gluon mode t —*■ eg. A similar study was performed for the FCNC 
top quark decays into Higgs bosons in a general two-Higgs-doublet model (2HDM) [TU] . 

The low observed rates of FCNC processes among the SM fields strongly suggest that 
their FCNC couplings must be suppressed at the tree-level in any extension of the SM, in 
particular in the 2HDM models. As is well known, there are two canonical strategies to get 
rid of these tree-level couplings in that kind of models jllj . In Type I 2HDM (also denoted 
2HDM I) one Higgs doublet, $1, does not couple to fermions at all and the other Higgs 
doublet, $2, couples to fermions in the same manner as in the SM. In contrast, in Type II 
2HDM (also denoted 2HDM II) one Higgs doublet, $1, couples to down quarks (but not to 
up quarks) while $2 does the other way around. Such a coupling pattern is automatically 
realized in the framework of supersymmetry (SUSY), in particular in the MSSM, but it 
can also be arranged in non-supersymmetric extensions if we impose a discrete symmetry, 
e.g. $1 — > — $1 and $2 — > +^2 (or vice versa) plus a suitable transformation for the 
right-handed quark fields. In an analysis of FCNC top quark decays in 2HDM extensions 
of the SM [Tmir2] it was proven that while the maximum rates for t eg were one order 
of magnitude more favorable in the MSSM than in the 2HDM, the corresponding rates 
for t — *> ch° were comparable both for the MSSM and the general 2HDM, namely up to 
the 10~ 4 level and should therefore be visible both at the LHC and the LC [TUllT2*i rEIj . 

Similarly, one may wonder whether the FCNC decays of the Higgs bosons themselves 
can be of some relevance. Obviously the situation with the SM Higgs is essentially hope- 
less, so again we have to move to physics beyond the SM. Some work on these decays, 
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performed in various contexts including the MSSM, shows that these effects can be impor- 
tant (T3JE3HII1HZ1IIE! • This could be expected, at least for heavy quarks in the MSSM, 
from the general SUSY study (including both strong and electroweak supersymmetric 
effects) of the FCNC vertices htc (h = h ,H° ,A°) made in Ref. |S]. However, other 
frameworks could perhaps be equally advantageous. Here we are particularly interested 
in the FCNC Higgs decay modes into top quark within a general 2HDM, which have not 
been studied anywhere in the literature to our knowledge. It means we restrict to Higgs 
bosons heavier than m t . From the above considerations, and most particularly on the 
basis of the detailed results obtained in [T01H2] one may expect that some of the decays 
of the Higgs bosons 

h^tc, h^ic (h = h°,H°,A°) (1) 

in a general 2HDM can be substantially enhanced and perhaps can be pushed up to 
the visible level, particularly for h° which is the lightest CP-even spinless state in these 
models [TT] . This possibility can be of great relevance on several grounds. On the 
one hand the severe degree of suppression of the FCNC Higgs decay in the SM obviously 
implies that any experimental sign of Higgs-like FCNC decay would be instant evidence 
of physics beyond the SM. On the other hand, the presence of an isolated top quark in 
the final state, unbalanced by any other heavy particle, is an unmistakable carrier of the 
FCNC signature. Finally, the study of the maximum FCNC rates for the top quark modes 
(JTJ) within the 2HDM, which is the simplest non-trivial extension of the SM, should serve 
as a fiducial result from which more complicated extensions of the SM can be referred 
to. Therefore, we believe there are founded reasons to perform a thorough study of the 
FCNC Higgs decays in minimal extensions of the Higgs sector of the SM and see whether 
they can be of any help to discover new physics. 

The paper is organized as follows. In Section 2 we summarize the 2HDM interactions 
most relevant for our study and estimate the expected FCNC rates of the Higgs decays 
in the SM and the general 2HDM. In Section 3 a detailed numerical analysis of the 
one-loop calculations of the FCNC decay widths and production rates of FCNC Higgs 
events is presented. Finally, in Section 4 we discuss the reach of our results and its 
phenomenological implications, and deliver our conclusions. 

2 Expected branching ratios in the SM and the 2HDM 

Before presenting the detailed numerical results of our calculation, it may be instructive 
to estimate the typical expected widths and branching ratios (BR) both for the SM 
decay H SM — > tc and the non-standard decays (JIJ in a general 2HDM. This should be 
especially useful in this kind of rare processes, which in the strict context of the SM are 
many orders of magnitude out of the accessible range. Therefore, one expects to be able 
to grossly reproduce the order of magnitude from simple physical considerations based 
on dimensional analysis, power counting, CKM matrix elements and dynamical features. 
By the same token it should be possible to guess at the potential enhancement factors in 
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the 2HDM extension of the SM. In fact, guided by the previous criteria the FCNC decay 
width of the SM Higgs of mass tub into top quark is expected to be of order 

T(H SM - f g) ~ (J^Y \V t tV bc fa 3 w m H (Af M ) 4 ~ f^j a w G 2 F m H mi, (2) 

where Gp is Fermi's constant and aw = g 2 /^, g being the SU{2) L weak gauge coupling. 
We have approximated the loop form factor by just a constant prefactor. Notice the 
presence of Xf ~ rrib/Mw, which is the SM Yukawa coupling of the bottom quark in 
units of g. The fourth power of Xf M in (J2J gives the non-trivial suppression factor remi- 
niscent of the GIM mechanism after summing over flavors. Since we are maximizing our 
estimation, a missing function related to kinematics and polarization sums, F{m t lmn) ~ 
(1— m 2 /m 2 H ) 2 , has been approximated to one. To obtain the (maximized!) branching ratio 
it suffices to divide the previous result by T(H SM — > bb) ~ aw (X§ AI ) mu ~ Gf mu m\ 
to obtain 

/ IT/ |\2 

T SM . -\ I \ y bc\ \ n 2 in-13 



BR{H bM - * c) ~ «w G F m fe 2 ~ 1(T 13 , (3) 

with VJ, C = 0.04, m& = 5 GeV. In general this 5i? will be even smaller, specially for higher 
Higgs boson masses (run > 2 Mw) for which the vector boson Higgs decay modes H SM — > 
W + W~(Z Z) can be kinematically available and become dominant. In this case it is easy 
to see that BR(H SM — > tc) will be suppressed by an additional factor of m 2 /m 2 H) which 
amounts at the very least to two additional orders of magnitude suppression, bringing it to 
a level of less than 10~ 15 . Already the optimized branching ratio (J5J) will remain invisible 
to all foreseeable accelerators in the future! To obtain the corresponding maximized 
estimation for the 2HDM we remind the following basic interaction Lagrangians within 
the general 2HDM. The charged Higgs interactions with fermions read 

C { £ tb = -JY^— H~b [m t cot (3P R + m b a s P L ] t + h.c. (4) 
v 2 Mw 

We use third-quark-family notation as a generic one; Pl,r = (1/2) (1 =F 7s) are the chiral 
projection operators on left- and right-handed fermions; j = I, II runs over Type I 
and Type II 2HDM's. A most relevant parameter here is tan/?, the ratio between the 
vacuum expectation values of the two Higgs doublet models [TT]. In the Lagrangian 
above, the parameter a,j is such that ai = — cot (3 and an = +tan/3. For the neutral 
Higgs interactions, the necessary pieces of the Lagrangian are the following: 



A3) _ ~9 m b 



'hqq 



2 Ml 



W 



h u i cosa I o | smn 



— sin a \ cos a 



sin (3 
cos (3 



b H 075 bA 

2M W 15 



+ / J \h° cosa + H sinal t+ J m J *y s tA° , (5) 
2 M w sin (3 1 J 2 M w tan (3 15 ' v ; 

where the upper row is for Type I models and the down row is for Type II. Here a is 
the mixing angle in the CP-even Higgs sector. Apart from the interactions with fermions, 
there is a set of potentially very relevant trilinear self-interactions of Higgs bosons in a 
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M w lm20 K m H± m 2 A0 + 2 m H°) Sm2 P COS (P <*) + 

lvlyv bill Zi LJ L v J--1 N 

+ {m 2 A0 — m 2 H0 ) cos2/?sin(/3 — a)] 




M w sm2/3 K m ff± m l« + 2^o)sm2/?sm(/? a) + 
+ (m^„ — cos 2/5 cos(/3 — a)] 


h°h°H° 


-S3 [(2m 2 h0 +m 2 H0 )sm2a- 
— m^ (3 sin 2a — sin 2/3)] 


A°A°H° 


2Afwsm2/3 K« sm 2/3 cos(/3 a)+ 
+2(m^ — ^^o) cos 2/? sin(/3 — a)] 


A°A°h° 


2M/sin2/3 K° sm2/?sm(/3 a)+ 
+2(m^ — m^ ) cos 2(3 cos(/3 — a)] 



Table 1: The needed trilinear self-couplings of the Higgs bosons in the 2HDM within the 
framework of Ref . jTUJ . 



general 2HDM. These are summarized in Tabled*. Notice that the trilinear couplings in 
Tableware common to the 2HDM I and 2HDM II, and they are quite different from the 
ones obtained in the MSSM, whose pure gauge structure is enforced by the underlying 
sup er sy mmet ry. 

Let us now first assume large tan (3 and restrict to Type II models. From the interaction 
Lagrangians above it is clear that we may replace Af M — > X^ M tan/3 in the previous 
formulae for the partial width. Moreover, the leading diagrams in the 2HDM contain the 
trilinear Higgs couplings A#+ h~ h- Therefore, the maximum BR associated to the FCNC 
decays Q in a general 2HDM II should be of order^ 

B R n (h^tc)^ (^^j a w G F m 2 b tzn 2 (3\ 2 H+H - h , (6) 

where A#+ h- k is defined here in units of g and dimensionless as compared to Table H 
Clearly a big enhancement factor tan 2 (3 appears, but this does not suffice. Fortunately, 
the trilinear couplings A#+ #- h for h = h°, H° (but not for h = A ) carry two additional 
sources of potential enhancement (Cf. Table which are absent in the MSSM case. Take 
e.g. h°, then we see that under appropriate conditions (for example, large tana and large 
tan/?) the trilinear coupling behaves as Xh+ h- h° ~ ( m h° ~ m A°) tan/3/(Mw m#±), and 
in this case 

BR n {h° - t c) ~ ( M") 2 aw Gp m l tan 4 (3 f ™f ~ ' . (7) 
\16n 2 J \ M w m H ± J 

So finally BR n (h° — > tc), and of course BR n (h° — > tc), can be augmented by a huge 
factor tan 4 (3 times the square of the relative splitting among the CP-even Higgs decaying 

* Tabled extends the one in Ref. [TOj, and therefore we also assume here A5 — \§ in the 2HDM Higgs 
potential [H]. We have omitted the couplings with the Goldstone bosons. Although the latter have 
been included in the calculation, their potential enhancement is far less significant than in the case of the 
Higgs trilinear couplings. 

'Here we have normalized the BR with respect to the h — > bb channel only, because the gauge boson 
modes will be suppressed in the relevant FCNC region, Cf. Section 3. 
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boson mass and the CP-odd Higgs mass. Since the neutral Higgs bosons do not participate 
in the loop form factors (Cf. Fig. 1 Ref. (101); it is clear that various scenarios can be 
envisaged where these mass splittings can be relevant. In the next section this behaviour 
will be borne out by explicit calculations showing that h — > tc can be raised to the visible 
level in the case of the Type II model. As for the Type I model the Higgs trilinear coupling 
enhancement is the same, but in the charged Higgs Yukawa coupling all quarks go with 
a factor cot /3; hence when considering the leading terms in the loops that contribute one 
sees that in the corresponding expression (JBJ the term m 2 tan 2 j3 is traded for m t m c cot 2 0, 
which is negligible at high tan /3. Both sources of enhancement are needed, and this feature 
is only tenable in the 2HDM II. Of course one could resort to the range tan /3 -C 1 for the 
Type I models, but this is not theoretically appealing. For example, for tan/5 < 0.1 the 
top quark Yukawa coupling g t = g m t /(\/2 M w sin/5), which is present in the interaction 
Lagrangians above, is pushed into the non-perturbative region g 2 /167r 2 > 1 and then 
our calculation would not be justified. And what is worse: for the 2HDM I we would 
actually need tan/3 < O(10~ 2 ) to get significant FCNC rates! In short, we consider that 
BR^h — > tc + ic) is essentially small (for all h), and that these decays remain always 
invisible to speak of. Hereafter we abandon the study of the decays (JTJ) for the 2HDM I 
and restrict ourselves to the general 2HDM II. 

3 Numerical analysis 

Let us now substantiate the previous claims and provide the precise numerical results of 
the full one-loop calculation of BR n {h — > tc + ic)* as well as of the LHC production 
rates of these FCNC events. We shall closely follow the notation and methods of Refs. 
[H1EI1E]- We refer the reader to these references for more details. See also ^Tj for 
basic definitions in the general 2HDM framework and [TTH l2"U] for calculational techniques 
and further illustration of the tan /3 enhancement in other relevant Higgs processes both 
in the MSSM and the 2HDM. In what follows we shall limit ourselves to present the 
final results of our numerical analysis together with a detailed discussion, interpretation 
and phenomenological application. We have performed the calculations with the help 
of the numeric and algebraic programs FeynArts, FormCalc and Loop Tools [21]. The 
calculation must obviously be finite without renormalization, and indeed the cancellation 
of UV divergences in the total amplitudes was verified explicitly. 
The input set for our numerical analysis is given by the data row 

(mho, mH0,TnA°,mH±, tana, tan /3) (8) 

made out of six independent parameters in the general 2HDM. Remaining inputs as 
in [22] ■ In practice there are some phenomenological restrictions on the data (jSJ) which 
were already described in |TU] and references therein, particularly [23]- It suffices to say 
here that the Higgs boson masses in (jHJ) are much less restricted than in the MSSM. 
Actually, the direct experimental searches put a soft bound of ~ 70 GeV for the lightest 

*Herc and throughout we use the notation BR H (h° -> tc + ic) = BR n (h° -> ic) + BR n (h° ic). 
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Figure 1: (a) BR n {h Q — > tc + ic) versus m^±; (b) Idem, versus m h o; (c) The production 
cross-section (in pb) of h° at the LHC versus its mass; (d) 5p mDM versus m^o, see the 
text. In these figures, when a parameter is not varied it is fixed as in eq.(JTTJ). 



neutral Higgs boson and around 80 GeV for the charged one H . These are much more 
relaxed than the bound of 114.1 GeV placed on the SM Higgs boson gU- In the MSSM, 
direct LEP 200 searches set a limit of > 90 GeV for m^o |24j and a similar limit for m^o 
which translates (approximately) to ~ 120 GeV due to the supersymmetric Higgs 
mass relations. Moreover, in contrast to the general 2HDM, in the MSSM case there 
is a theoretical upper bound m/jo < 122, 135 GeV (for minimal and maximal top-squark 
mixing, respectively) |25j . Concerning tan/3, in a general 2HDM it is restricted only, by 
perturbativity arguments, to lie within the approximate range 

0.1 < tan/? < 60. (9) 

Strictly speaking the upper bound could be larger, but as in [TU] we shall stick to the 
above moderate range. In practice, since Type I models are not considered, the effective 
range for our calculation will be the high tan (5 end of (jHJ). For comparison, in the MSSM 
case there are additional phenomenological restrictions that bring the lower bound on 
tan/3 to roughly 2.4 for the so-called maximal m° h scenario, and 10.5 for the so-called no 
mixing scenario [21], the upper bound being the same. Furthermore, in contrast to the 
MSSM, the parameter tan a is free in the 2HDM. 
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Figure 2: Contour lines in the (m/j±, m/ l o)-plane for the branching ratios (2HDM II case) 
(a) BR H (h° -f t c + t c) and (b) BR H (H° ^tc + ic) assuming 5p mDM at 1 a. 

Apart from the direct limit uih± > 80 GeV from LEP 200, there are more stringent 
restrictions on the charged Higgs boson mass. Although they are indirect and in this sense 
not so indisputable, it is wise to assess their impact in our analysis. These restrictions 
emerge from the virtual contributions to b — > sj that have been computed at the NLO 
in QCD [26J. From these calculations and from the experimental limits on the radiative 
B-meson decays, these authors can place a lower bound on the 2HDM II charged Higgs 
mass of m H ± > 350 GeV We will, in general, apply this lower bound, but since the 
direct and indirect searches complement each other, we may disregard it in some parts 
of our numerical analysis. The parts where we deviate will be clear in the text. We 
recall that at the moment the MSSM may escape this indirect bound because the positive 
charged Higgs virtual contributions to b — > 37 can be compensated for by negative stop 
and chargino loops, if they are not too heavy. Therefore, in the MSSM the charged Higgs 
can stay relatively light, > 120 GeV, just to comply with the aforementioned LEP 
200 bounds on m A o [21 . 

Finally, we will restrict the input data row (JSJ) with the important constraint from Sp, 
extensively used in |H1GII] for the FCNC top quark decays in the MSSM and the 2HDM. 
Namely, the one-loop corrections to the p-parameter (i.e. the ratio between the charged 
and neutral current Fermi constants) are bound within |5p 2HDM | < 0.1%. To be precise, 
the latter is the extra effect that 5p can accommodate at one standard deviation (ler) 
from the 2HDM fields beyond the SM contribution |H]. This is a stringent restriction 
that affects the possible mass splittings among the Higgs fields of the 2HDM, and its 
implementation in our codes does severely prevent the possibility from playing with the 
Higgs boson masses to artificially enhance the FCNC contributions. 

With these restrictions in mind we have computed the number of FCNC Higgs decay 
events into top quark at the LHC: 



pp^ h + X -*tc(tc) +X (h = h°,H°,A°) . 



(10) 
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Figure 3: (a) Contour lines in the (m#±, m h o)-pla.ne for the maximum number of light 
CP-even Higgs FCNC events h° —>■ tc+ic produced at the LHC for 100 fb~ l of integrated 
luminosity; (b) Contour lines showing the value of m^o that maximizes the number of 
events; (c) As in (a) but within 5p 2UDM at 3 a; (d) As in (b) but with 5p 2UDM at 3 a. 



The necessary cross-sections to compute the production of neutral Higgs bosons at this 
collider, including all known QCD corrections, have been computed by adapting the codes 
HIGLU 1.0 and HQQ 1.0 j23 - originally written for the MSSM case [H] to the general 
case of the 2HDM§. Folding the cross-sections with the one- loop branching ratios of the 
processes (JTJ) we have obtained the number of FCNC Higgs decay events at the LHC. Let 
us first consider the branching ratios themselves. In Fig.[TJi,b we show BR n (h° — > t c+i c) 
for the lightest CP-even state. In particular, Fig. shows BR H (h° — ► tc + ic) versus 
the charged Higgs mass mn±. In this figure we fix the values of the parameters in (JHJ) 

§We have used the default parton distribution functions and renormalization/factorization scales used 
in these programs, namely GRV94 with = /if = rrih for HIGLU, and CTEQ4L with = /if = •*/§ = 

^{Ph+PQ+PQ? for HQQ. 



9 




Figure 4: Contour lines a/ir = const, (a is the mixing angle in the CP-even sector) 
corresponding to Fig. El for 5p mDM at (a) 1 a and (b) 3 a . 



which are not varying as follows: 

(m h o = 350 GeV, m H o = 600 GeV, m A o = 550 GeV, m H ± = 375 GeV, ( . 

tana = 30, tan/? = 60) 

After crossing a local maximum (associated to a pseudo-threshold of the one-loop vertex 
function involving the h° H + H~ coupling) the subsequently falling behavior of the BR 
with mn± clearly shows that the previously discussed bounds on m#± are quite relevant. 
The branching ratio, however, stays within 10~ 6 — 10~ 5 for a wide range of heavy charged 
Higgs masses extending up to m H ± < 600 GeV in Fig. [TJi. Hence, for m H ± heavy enough 
to satisfy the indirect bounds from radiative 5-meson decays the maximum BR is 
still sizeable. In Fig.^i the production rate of h° bosons at the LHC is shown as a function 
of m/jo, for fixed parameters (jHJ)- The production cross-sections for the subprocesses 

ff 3 _> h° + X, gg,qq-+ h° + QQ, (12) 

contributing to (|10|) in the case of the light CP-even Higgs h° are explicitly separated in 
Fig. Ht. The gluon-gluon fusion process proceeds at one-loop and the h°QQ associated 
production proceeds at tree- level [29^ Similar subprocesses and results apply for H° 
and A production. At large tan /3 and the larger the Higgs boson masses the particular 
associated production mechanism with the bottom quark, Q = b, i.e. h° bb, becomes 
dominant by far. All other mechanisms for Higgs boson production in Type II mod- 
els |2l^l2l?|,l3Ti] , like vector-boson fusion (which contributes also to h°QQ when Q are light 
quarks), vector-boson bremsstrahlung (qq — > hV) and associated tt production, are sub- 
dominant at large tan (3 and can be neglected for our purposes. Admittedly, some of these 
mechanisms can be relevant for Higgs boson production in the case of the Type I 2HDM 
at low tan /3, but we have already warned that the corresponding FCNC branching ratios 
are never sufficiently high. 
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Figure 5: (a) Contour lines in the (m#±, m#o)-plane for the maximum number of heavy 
CP-even Higgs FCNC events H° -> tc + ic (2HDM II case) produced at the LHC for 
100 fb^ 1 of integrated luminosity; (b) Corresponding contour lines for m^o; (c) As in (a) 
but within 5p mDM at 3 a; (d) As in (b) but with 5p 2UDM at 3 a. 



The control over 5p 2HDM is displayed in Fig. ^i. Recall that 5p is not sensitive to 
the mass splitting between m h o and m H o, because of CP- conservation in the gauge boson 
sector, but it does feel all the other mass splittings among Higgs bosons, charged and 
neutral. A more systematic search of BR values in the parameter space is presented in 
Figs. |2K ; b corresponding to BR n (h° —>■ tc + ic) and BR n (H° — > t c + t c) respectively. 
Here we have scanned independently on the parameters (jHJ) while holding the 5p 2HDM 
bound at 1 a. The contour lines in these figures represent the locus of points in the 
(mji±, m/jo)-plane giving maximized values of the BR in the 2HDM II. Let us remark 
that the highest value of tan/3 is always preferred, and therefore all these contour lines 
correspond to tan (3 = 60. 

In practice, to better assess the possibility of detection at the LHC, one has to study the 



11 




Figure 6: Contour lines a/n = const, as in FiglH but for the heavy CP-even Higgs. 

production rates of the FCNC events. These are determined by combining the production 
cross-sections of neutral 2HDM II Higgs bosons at the LHC and the FCNC branching 
ratios. If we just adopt the mild LEP bound m H ± > 80 GeV and let m H ± approach the 
maximum in Fig. then the BR can be as large as 10~ 3 and the number of FCNC 
events can be huge, at the level of ten thousand per 100 fb -1 of integrated luminosity. 
But of course the region near the maximum is too special. Moreover, if we switch on the 
above mentioned indirect bound from 6 — > s 7 , then the typical BR is much smaller 
(of order 10~ 5 ) and the number of events is reduced dramatically, at a level of hundred 
or less for the same integrated luminosity. On the other hand it may well happen that 
there are regions of parameter space where BR ~ 10 -5 (see Fig. EJ) but the production 
cross-section is too small because the decaying Higgs boson is too heavy. Therefore, it is 
the product of the two quantities that matters. 

The systematic search of the regions of parameter space with the maximum number 
of FCNC events for the light CP-even Higgs is presented in the form of contour lines in 
the multiple Fig. El For instance, each isoline in Figs. |3]a,c corresponds to a fixed number 
of produced FCNC events at the LHC while keeping the value of 5p mDM within 1 a or 3 a 
respectively of its central experimental value. When scanning over the parameter space 
(JBJ) we have found again that tan (3 is preferred at the highest allowed value (tan (3 = 60) - 
for Type II models. We have also determined (see Figs.|3]b,d) the corresponding contour 
lines for m^o associated to these events. The m^o-lines are important because the FCNC 
processes under consideration are sensitive to the mass splittings between m^o and the 
corresponding decaying Higgs boson, see e.g. eq.(J7J) and Tabled] The combined figures 
|3]a-d are very useful because they give a panoramic view of the origin of our results in 
the parameter space. To complete the map of the numerical analysis we provide Fig. 0] 
in which we have projected the contour lines of the CP-even mixing angle a associated 
to the previous plots. For a given contour line a/n = const., the set of inner points have 
a value of a /it smaller than the one defined by the line itself. In particular, the large 
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domains in Figs. Hk,b without contour lines correspond to a/n > 0.4 and so to relatively 
large (and positive) tana. There are a few and small neighborhoods where the FCNC 
rates for h° can be sizeable also for small tana. 

Knowing that high tan a is generally preferred by h — > tc + ic, and noting from 
Fig-El that large mass splittings between m^o and m^o are allowed, we find that the 
trilinear coupling X H + h- h° can take the form X H + h- h° ~ ( m ho — m \o ) tan f3 / [Mw m H± ) • 
Hence it provides a substantial additional enhancement beyond the tan/3 factor. One 
can check from the approximate formula (J7J) that the maximum FCNC branching ratios 
BR n (h° — > tc+ic) can eventually reach the 10~ 5 level even in regions where the charged 
Higgs boson mass preserve the stringent indirect bounds from radiative B-meson decays 
|26j . These expectations are well in agreement with the exact numerical analysis presented 
in Fig. |21 thus showing that eq. (J7|) provides a reasonable estimate, and therefore a plausible 
explanation for the origin of the maximum contributions. As a matter of fact, we have 
checked that the single (finite) Feynman diagram giving rise to the estimation ((7j) - the 
one-loop vertex Feynman diagram with a couple of charged Higgs bosons and a bottom 
quark in the loop - reproduces the full result with an accuracy better than 10% for 
tan/? > 10 — 20. At lower tan/3 values large deviations are possible but, as warned 
before, eq. (J7|) is expected to be valid only at large tan/3. Furthermore, for low values of 
tan/3 < 20 the FCNC BRs are too small to be of any phenomenological interest. The 
exact numerical analysis is of course based on the full expression for the branching ratio 

R/?W^f,xU- T(h° ^ tc + ic) 

{ 4 } Y(h^bb) + Y(h^tt) + Y(h^VV) + Y(h^HH)' K } 

where all decay widths in the denominator of this formula have been computed at the 
tree-level in the 2HDM II, since this provides a consistent description of eq. (|T^|) at leading 
order. Here we have defined 

Y(h° -+VV) = T(h° -> W + W~) + Y(H° -+ZZ), (14) 



Y(h° -+HH) = T(h° -> A A ) + Y(h° -> H + H~) . (15) 

We disregard the loop induced decay channels, since they have branching ratios below the 
percent level all over the parameter space. The r-lepton decay channel is also neglected, 
since it is suppressed by a factor of 0(1O~ 2 ) with respect the 66-channel in the whole 2HDM 
parameter space. In general the effect of the gauge boson channels h° —>■ W + W~~ , ZZ 
in the BR fllMjl is not so important as in the SM, actually for (3 = a they vanish in the 
h° case because they are proportional to sin 2 (/3 — a). This is approximately the case 
for large tana and large tan/3, the dominant FCNC region for h° decay (Cf. Fig0Ji 
and Eh). In this region, the mode h° — > tt is, when kinematically allowed, suppressed: 
BR(h° -> tt) oc cos 2 a I sin 2 /3 -> (Cf. Eq. ©). On the other hand there are domains in 
our plots where the decays h° — > H + H~ and h° — > A A are kinematically possible and 
non- (dynamically) suppressed. Indeed, this can be checked from the explicit structure 
of the trilinear couplings h° H + H~ and h° A A in Table Q in the dominant region for 
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the decays h° — > tc + tc both of these couplings are tan /5-enhanced. Nevertheless the 
decay h° — > A A is only possible for m^o < m h o /2 , and since the optimal FCNC 
regions demand the largest possible values of m A o, this decay is kinematically blocked 
there. On the other hand the mode h — > H + H is of course allowed if we just take the 
aforementioned direct limits on the 2HDM Higgs boson masses. But it is never available 
if we apply the indirect bound from b — > s 7 on the charged Higgs mass mentioned above, 
unless m h o > 2m H ± > 700 GeV , in which case h° is so heavy that its production cross- 
section is too small for FCNC studies to be further pursued. 

The corresponding results for the heavy CP-even Higgs boson are displayed in Figs. 
and El The exact formula for the BR in this case reads 

BR n (H° ^tc + ic) = T(H°-+tc + ic) 

1 ' T(H° ^bb)+T(H° ^ti) + T(H° ^VV) + T(H° -> H H) ' 

(16) 

where we have defined 

T(H° -»• VV) = T(H° -> W + W~) + T(H° -»• Z Z) , (17) 



T{H° ^HH)= T(H° -> h° h°) + T(H° -> A A ) + T(H° -> F") . (18) 

From the contour lines in Fig.|S^,c it is patent that the number of FCNC top quark events 
stemming from H° decays is comparable to the case of the lightest Higgs boson. However, 
Fig. Et,b clearly reveals that these events are localized in regions of the parameter space 
generally different from the h° case, namely they prefer tana ~ 0. Even so, there are some 
"islands" of events at large tana. This situation is complementary to the one observed 
for h° in Fig.0] However, in both cases these isolated regions are mainly concentrated in 
the segment < 350 GeV. Therefore, if the bound on m#± from b — >• s 7 is strictly 

preserved, it is difficult to find regions of parameter space where the two CP-even states 
of a general 2HDM II may both undergo a FCNC decay of the type (0). 

In the dominant regions of the FCNC mode H° — > tc (where tan a is small and 
tan/? is large), the decay of H° into the tt final state is suppressed: BR(H — ► tt) oc 
sin 2 a/ sin 2 /3 — * 0. In the same regions the gauge boson channels in (|16|) are suppressed 
too because T(H° — > W + W~ , ZZ) oc cos 2 (/? — a). In principle the heavy CP-even Higgs 
boson H° also could (as h°) decay into A A and H + H . But there is a novelty here 
with respect to the h° decays, in that there could be regions where H° could decay into 
the final state h° h°. This contingency has been included explicitly in eq. ()18|) . However, in 
practice, neither one of these three last channels is relevant in the optimal FCNC domains 
of parameter space. First, the decay H° — > h° h°, although it is kinematically possible, 
is dynamically suppressed in the main FCNC region for H°. This can be seen from 
Table ^ where the trilinear coupling H° h° h° becomes vanishingly small at large tan f3 
and small tana. Second, the coupling H°A°A° in Table[T]is non-suppressed in the present 
region, but again the mode H° — > A A is kinematically forbidden in the optimal FCNC 
domains because the latter favor large values of the CP-odd mass (see Fig.GJ>,d). Third, 
although in these domains the decay H° — > H + H~ is also non-dynamically suppressed 
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(see the corresponding trilinear coupling in Table [TJ, it becomes kinematically shifted to 
the high mass range m H o > 700 GeV if we switch on the indirect bound from b — > 57. 
Obviously, in this latter case the H° production cross section becomes too small and the 
FCNC study has no interest. All in all the contributions from (|I3 j) .(|I5 j) .([r7 j) and (fTHjl are 
irrelevant for m^o, m#o < 700 GeV as their numerical impact on B H (h ,H° — > tc + tc) 
is negligible. Our formulae (fT3j) and (fTtjj) do contain all the decay channels and we have 
verified explicitly these features. 

As remarked before, in general the most favorable regions of parameter space for the 
FCNC decays of h° and H° do not overlap much. The trilinear Higgs boson self-couplings 
in Table ^ (also the fermionic ones) are interchanged when performing the simultaneous 
substitutions a — > n/2 — a and m/jO — > m#o [TU]. Furthermore, the LHC production rates 
of the neutral Higgs bosons fall quite fast with the masses of these particles, as seen e.g. 
in Fig. for the h° state. As a consequence that exchange symmetry on the branching 
ratios does not go over to the final event rates, so in practice the number of FCNC events 
from H° decays are smaller (for the same values of the other parameters) as compared 
to those for h°; thus H° requires e.g. lighter charged Higgs masses to achieve the same 
number of FCNC events as h°. As for the CP-odd state A , we have seen that it plays an 
important indirect dynamical role on the other decays through the trilinear couplings in 
Table [TJ but its own FCNC decay rates never get a sufficient degree of enhancement due 
to the absence of the relevant trilinear couplings, so we may discard it from our analysis. 

We notice that this picture is consistent with the decoupling limit in the 2HDM: 
for a — > (3, the heaviest CP-even Higgs boson (H°) behaves as the SM Higgs boson, 
whereas h° decouples from the electroweak gauge bosons and may develop enhanced 
couplings to up and down-like quarks, depending on whether tan/3 is small or large 
respectively; in the opposite limit (a — ► (3 — vr/2), it is h° that behaves as H , while H° 
decouples from gauge bosons and may develop the same enhanced couplings to quarks 
as h° did in the previous case. Indeed these are the situations that we find concerning 
the FCNC decay rates. We recall that the numerical results presented in our figures 
correspond to an integrated luminosity of 100 fb~ l . However, the combined ATLAS and 
CMS detectors might eventually accumulate a few hundred inverse femtobarn pHl l32|. 
Therefore, hopefully, a few hundred FCNC events could eventually be collected in 
the most optimistic scenario. Actually, the extreme rareness of these events in the SM 
suggests that if only a few of them could be clearly disentangled, it should suffice to claim 
physics beyond the SM. 

4 Discussion and conclusions 

Detection strategies at the CERN-LHC collider for the search of the SM Higgs boson, 
and also for the three spinless fields of the MSSM Higgs sector, have been described 
in great detail in many places of the literature jTTJEHE21E3Ell, but not so well for 
the corresponding charged and neutral Higgs bosons of the general 2HDM. The result 
is that the discovery of the SM Higgs boson is guaranteed at the LHC in the whole 
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presumed range 100 GeV < tuh < 1 TeV. However, the discovery channels are different 
in each kinematical region and sometimes the most obvious ones are rendered useless. For 
example, due to the huge irreducible QCD background from bb dijets, the decay mode 
H SM — s> bb is difficult and one has to complement the search with many other channels, 
particularly H SM — > 77 pTl) l32| . We have shown in this work that there are scenarios in 
the 2HDM parameter space where alternative decays, like the FCNC modes h° — > tc + ic 
and H° -> t c + i c, can also be useful. For instance, in the h° case, this situation occurs 
when tan/3 and tana are both large and the CP-odd state is much heavier than the 
CP-even ones. The potential enhancement is then spectacular and it may reach up to 
ten billion times the SM value BR(H SM — > tc) ~ 10~ 15 , thereby bringing the maximum 
value of the FCNC branching ratio BR(h° —>■ tc) to the level of ~ 10 -5 . As a matter of 
fact, the enhancement would be much larger were it not because we eventually apply the 
severe (indirect) lower bound on the charged Higgs mass from b — > sj Although 
these decays have maximal ratios below BR(h —>■ 77) ~ 10~ 3 , they should be essentially 
free of QCD background 

While in the MSSM almost the full (m^o, tan/?)-parameter space is covered, with 
better efficiency at high tan /3 though, we should insist that within the general 2HDM the 
tagging strategies are not so well studied and one would like to have further information 
to disentangle the MSSM scenarios from the 2HDM ones. Here again the study of the 
FCNC Higgs decays can play a role. Of course the statistics for the FCNC Higgs decays 
is poor due to the weakness of the couplings and the large masses of the Higgs bosons 
to be produced. However, in the favorable regions, which are generally characterized 
by large values of tan/3 and of tan a, one may collect a few hundred events of the type 
(Q)- mainly from h° - in the high luminosity phase of the LHC. As we have said, this 
is basically due to the enormous enhancement that may undergo the FCNC decay rates, 
but also because in the same regions of parameter space where the BR's are enhanced, 
also the LHC production rates of the Higgs bosons can be significantly larger (one order 
of magnitude) in the 2HDM II as compared to the SM. 

Interestingly enough, in many cases one can easily distinguish whether the enhanced 
FCNC events stem from the dynamics of a general, unrestricted, 2HDM model, or 
rather from some supersymmetric mechanisms within the MSSM. This is already obvious 
from the fact that the ranges of neutral and charged Higgs boson masses in the 2HDM 
case can be totally incompatible with the corresponding ones in the MSSM. But there 
are many other ways to discriminate these rare events. For instance, in the 2HDM case 
the CP-odd modes A — > tc + ic are completely hopeless whereas in the MSSM they 
can be enhanced Using this information in combination with the masses 

of potentially detected Higgs bosons could be extremely useful to pinpoint the supersym- 
metric or non-supersymmetric nature of them. We may describe a few specific strategies. 
As it was first shown in Ref. [S], the leading SUSY-FCNC effects associated to the htc 
vertices (h = h° , H° , A ) come from the FCNC gluino interactions which are induced 

^Misidentification of 6-quarks as c-quarks in tb production might be a source of background to our 
FCNC events. However, to rate the actual impact of that misidentification one would need a dedicated 
simulation of the signal versus background, which is beyond the scope of this paper. 
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by potentially large misalignments of the quark and squark mass matrices These 
effects are not particularly sensitive to tan (3 and they can be very sizeable for both high 
and moderately low values of this parameter. This sole fact can be another distinguishing 
feature between FCNC events (QJ) of MSSM or 2HDM origin. If, for example, a few of 
these events were observed and at the same time the best MSSM fits to the electroweak 
precision data would favor moderate values of tan/3, say in the range 10 — 20, then it is 
clear that those events could originate in the FCNC gluino interactions but in no way 
within the context of the general 2HDM. In this respect it should be mentioned that the 
FCNC gluino couplings recently became more restricted from the low-energy meson data 
and will presumably become further restricted in the near future. The reason being 
that the same couplings are related, via SU(2) gauge invariance and CKM rotation, to 
those affecting the down-like quark sector, which will most likely become constrained by 
the increasingly more precise low-energy meson physics [HHtlHT] . In that circumstance the 
only source of FCNC Higgs decays in the MSSM will stem purely from the electroweak 
interactions within the super-CKM basis. Then, in the absence of these SUSY-QCD 
FCNC effects, we could judiciously conclude from the work of Ref. - in which both the 
SUSY-QCD and the SUSY electroweak contributions were computed for the htc vertices 
- that the FCNC rates in the MSSM should diminish dramatically (two to three orders 
of magnitude). In such case we can imagine the following "provocative" scenario. Sup- 
pose that the LHC finds a light neutral Higgs boson of mass < 140 GeV (suggestively 
enough, in a mass range near the MSSM upper bound for m h o\) and subsequently, or 
about simultaneously, a charged Higgs boson and another neutral Higgs boson both with 
masses around 400 GeV or more. At this point one could naively suspect that a MSSM 
picture out these findings is getting somehow confirmed. If, however, later on a few FCNC 
events (0) are reported and potentially ascribed to the previously discovered heavy neu- 
tral Higgs boson (presumably H°), then the overall situation could not correspond at all 
to the MSSM, while it could be perfectly compatible with the 2HDM II. Alternatively, 
suppose that the FCNC gluino couplings were not yet sufficiently restricted, but (still 
following the remaining hypotheses of the previous example) a third neutral Higgs boson 
(presumably A ) is found, also accompanied with a few FCNC events. Then this situation 
would be incompatible with the 2HDM II, and in actual fact it would put forward strong 
(indirect) evidence of the MSSM!! 

We should also mention that there are other FCNC Higgs decay modes, as for example 
h — > bs + bs, which could be, in principle, competitive with the top quark modes (JI)i. 
In some cases these bottom modes can be highly enhanced in the MSSM case [To^ ll6j. 
Actually, a more complete assessment of the FCNC bottom modes in the MSSM case - 
namely one which takes also into account the supersymmetric contributions to the highly 
restrictive radiative B-meson decays - shows that they are eventually rendered at a similar 
level of the top modes under study in most of the parameter space [T5j . 

To summarize, the FCNC decays of the Higgs bosons into top quark final states can 
be a helpful complementary strategy to search for signals of physics beyond the SM in the 
LHC. Our comprehensive numerical analysis shows that the FCNC studies are feasible for 
CP-even Higgs masses up to about 500 GeV. While the statistics of these FCNC decays 
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is of course poor, the advantage is that a few tagged and well discriminated events of this 
sort could not be attributed by any means to the SM, and therefore should call for various 
kinds of new physics. In this paper we have shown that a general 2HDM II is potentially 
competitive to be ultimately responsible for these FCNC decays, if they are ever found, 
and we have exemplified how to discriminate this possibility from the more restricted one 
associated to the MSSM. 
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